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We discuss the properties of superintegrable Hamiltonian systems, in particular those that admit 
separation of variables in cartesian coordinates. We show that the superintegrability of such poten- 
tials is equivalent to the isochronicity of the separated potentials. We use this fact to get a new 
insight into an old question about the relation between quantum and classical harmonic behavior. 
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H = pl/2 + p 2 y /2 + V a (x) + V b (y) 



(1) 



admit two independent second-order integrals, H a — 
pl/2 + V a (x) and H b = p%/2 + V b {y). They are therefore 
integrable. Those that admit an additional nontrivial 
integral of motion are maximally superintegrable since 
they have 2n — 1 integrals of motion, n being the number 
of degrees of freedom. Such systems were investigated in 

Ml. 

Theorem 1X1 relates the existence of this additional in- 
tegral of motion to the isochronicity of the separated 
Hamiltonians. It provides a new characterization of 
isochronous potentials in classical mechanics, which can 
easily be adapted to many more general contexts. We 
use this to propose a new definition of quantum harmonic 
behavior, based on quantum superintegrability. We then 
give examples that show that potentials that are har- 
monic according to this definition do not necessarily have 
equidistant spectra, which was the usual definition for 
quantum harmonic behavior, but still exhibit regular be- 
havior. The spectra of our examples can indeed be gener- 
ated from a finite number of states by the application of 
a creation operator. Therefore their emission spectrum 
is highly degenerate, which makes them physically very 
similar to the usual harmonic oscillator. 

Remark 1. We will assume through this paper that all 
classical potentials are C2 over the considered interval of 
the real axis. The results might be generalizable to C\. 

Definition 1. A Hamiltonian H a (x,p x ) and, by exten- 
sion, the potential V a (x) is isochronous if V a has a lo- 
cal minimum and if all its bounded trajectories have the 
same period. 

Remark 2. Since the potentials considered here are dif- 
ferentiable, and the existence of a local maximum would 
imply trajectories with infinite period, an isochronous 
potential can have only one local minimum, and the po- 
tential must go to infinity on the boundaries of the con- 
sidered domain. 

Theorem A. A Hamiltonian of the form Q with a local 
minimum is maximally superintegrable if and only if V a 
and V b are isochronous. 



Proof. The first part of the proof is straightforward. We 
know, from a theorem due to Nekhoroshev p£], that 
the bounded trajectories of a maximally superintegrable 
Hamiltonian are all closed. Since the motion along the x 
and y axis can be decoupled, we can consider separately 
the periods along each direction. If these periods varied, 
they would do it continuously with respect to initial con- 
ditions that are independent (i.e. the 2i(0) and Pi(0)). 
We could therefore choose initial conditions such that the 
ratio between the periods is irrational. In this case the 
two-dimensional motion would not be periodic, leading 
to a contradiction. 

We now have to demonstrate that the converse is true, 
i.e. that if V a and V b are isochronous the whole system 
is superintegrable. The idea of the proof is to trans- 
form each one-dimensional system into a harmonic os- 
cillator, and then use the two-dimensional integrals of 
the anisotropic harmonic oscillator to find, via the in- 
verse transformation, integrals of the initial system. We 
will find directly the appropriate transformation using 
the trajectories of the system. 

We can specify a point in the phase space of a one 
dimensional isochronous potential that has its minimum 
at x = by specifying the trajectory on which it lies, 
and the time it took for a particle moving along that 
trajectory to get there starting from the turning point 
with x > 0. Since there is only one such turning point, 
and since the system is isochronous, this transformation 
is well defined and one-to-one on the interval t E [0, T) , 
T being the period of the system. 

We therefore consider the transformation (x,p) — > 
(r, t), where r is the value of the turning point with x > 
and t the time for a particle to get to (x,p) from that 
turning point. We will now introduce a lemma that will 
be demonstrated below. 

Lemma 1. The transformation (x,p) — » (r, t) is a dif- 
feomorphism for r =/= 0. 

We can compose this change of variables with 



X = r cos t , P — r sin t , 



(2) 



which is a transformation from polar to cartesian coordi- 
nates and is continuous and differentiable, except at the 
origin. The transformation (x,p) — > (X,P) is therefore 
a diffeomorphism, except possibly at the origin. 
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The trajectories in the (X, P) plane have constant ve- 
locity, period T and are along circles. Therefore X and 
P obey the equations of a harmonic oscillator, 



X = P 



P=- 



X 



(3) 



We use this method for both {xi,pi) and {x2,P2), 
and since the periods in the two directions are com- 
mensurable, the resulting system expressed in terms of 
{Xl, X2, Pi, P2} exhibits the same motion, and therefore 
the same integrals as the anisotropic harmonic oscilla- 
tor with rational ratio. This is known || to admit three 
independent integrals of motion, {Qi,Q2,Qs}- 

Since the and Pj are differentiable functions of the 
Xi and The three Qi, expressed in terms of the {xi,pt} 
are independent integrals of motion of the initial Hamil- 
tonian. 

□ 

Let us now prove lemma ^ 

Proof. Differentiability of (x,p) with respect to time is 
given by the Hamilton equations, x = p and p = — V (x), 
The existence of the two derivatives with respect to r 
is given by the theorem of differentiable dependance on 
initial conditions (see e.g. jfj), applied to the Hamilton 
equations. 

In order to show that (x,p) — > (r,t) is also differen- 
tiable, we have to show that the determinant of the Ja- 
cobian matrix 



J = 



dx 


dx 












dp 




dt 


dr 





P(t,r) § 
-V'(x(t,r)) | 



does not vanish on a solution (x(t,r),p(t,r)) . Let us 
write 

x(t,r + v)=x v , p(t,r + v)=p v , x = x, p = p. 

The Jacobian reads 

dp Ox 
p— + V'(x)— = ]hap(p v — p)/v + V (x)(x v - x)/v 
Or or u->0 

= lim (p 2 v /2 - p 2 /2 + V'{x) (x v ~ x)) /v 
-\iMp v -pf/{2v) 

Since V'(x)(x v — x) ~ V(x v ) — V(x) and (p v —p) 2 /v ~ 
(dp/dr) 2 v ~ we have 

J = ]im(p 2 v /2 + V(x v )- (p 2 /2 + V(x)))/v 

v—>0 

= lim(F(r) - V(r + v))/v = -V'(r). 

This is nonzero everywhere except at r = 0. The proof 
of the lemma is complete. □ 

We now introduce a concept which will allow us to dis- 
cuss one-dimensional potentials without always referring 
to the two-dimensional Hamiltonian. 



Definition 2. A Hamiltonian H a = p 2 /2 + V a (x) is 
called 2£>-superintegrable if there exists a = p 2 /2 + 
Vb(y) such that H a + Ht is superintegrable. We shall 
call Hi,, Hamiltonian associated to H a . It is also 2D- 
superintegrable. 

Remark 3. If H a is isochronous, theorem 1X1 guarantees 
that it is 2£>-superintegrable. The fact that H a is 2D- 
superintegrable gives no guarantee that it is isochronous, 
though, unless both H a and its associated Hamiltonian 
have a local minimum. 

2_D-superintegrability is therefore a generalization of 
isochronicity since it includes potentials such as free mo- 
tion and the repulsive harmonic oscillator. It is therefore 
tempting to identify as "harmonic" 2D-superintegrable 
potentials in classical mechanics. 

Even though some questions regarding the indepen- 
dence of quantum integrals are not yet fully understood 
0, @ , the concept of superintegrability can be adapted 
to quantum mechanics in a straightforward manner, and 
this had been done since 0]. The adaptation is not so 
easy for isochronicity. Since well-known examples (such 
as the harmonic oscillator) have equidistant spectra, it 
was natural to assume that this prop erty should charac- 
terize quantum harmonic potentials [I(J, fill fl2| . It was 
shown in 0, ^| that the quantum equivalent to a clas- 
sical isochronous potential need not be equidistant. We 
will see that this is still the case when we ask that the 
quantum equivalent be 2£>-superintegrable. The quan- 
tum potential will nevertheless exhibit properties that 
make it similar to the harmonic oscillator, e.g. regarding 
the infinite degeneracy of the emission spectrum. There- 
fore it is tempting to identify quantum harmonic behav- 
ior and 2Z)-superintegrability. Studying the properties 
of 2D-superintegrable potentials will then help us in the 
search for a more "property-oriented" definition for quan- 
tum harmonic behavior. 

The formalism of creation-annihilation operators, de- 
fined here as operators A such that [H, A] — XA with A 
a constant, is useful in dealing with 2D-superintegrable 
potentials. Let us write H = H a {xi,pi) + H(,(x2,P2), 
and let Q{xi,X2,Pi,P2) be the additional integral of 
motion. If the commutator [H a ,Q] is nonzero, it is a 
new integral of H. This is likely to be a consequence 
of the definition of the independence of quantum inte- 
grals since it is the case in classical mechanics. There- 
fore the operator M : Q — > [H a , Q] is a nontrivial and 
nonzero linear operator on the vector space of the in- 
tegrals of motion of H. Therefore each nonzero eigen- 
value of M indicates the existence of a creation or an- 
nihilation operator for H a . Notice that the existence of 
creation-annihilation operators does not guarantee that 
the spectrum is equidistant. Consider for example the 



Hamiltonian H = p /2 



b/x 2 , which is classical 



and quantum 2_D-superintegrable for all values of a and 
b. Its spectrum and eigenf unctions are known In the 
quantum case we can find a creation-annihilation pair by 
the method just described, with |A| = 2\f2ah. 
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For b > 3/8ft 2 , the spectrum of this potential is equidis- 
tant. For — ft 2 /8 < b < 3ft 2 /8, though, the spectrum is 
given by 



where f(x) 



E = V2a(2k + 1 ± i>)h, 



where v = y/l + 86/ft 2 /2. Therefore, apart from the 
special case 6 = 0, the potentials with — l/8ft 2 < b < 
3/8ft 2 , have a spectrum organized in pairs separated by 
2\[2avh. Each pair is separated by 2v / 2aft- The creation- 
annihilation operators therefore skip a level every time 
they are applied. Finally notice that b is proportional 
to ft 2 , and therefore v does not depend on ft. Hence the 
energy levels are still grouped in pairs as ft — > 0. 

The method of dressing chains (see 0, Q|) can be 
used to construct Hamiltonians for which we can calcu- 
late explicitly the spectrum of energy. The idea is to con- 
sider a chain of one-dimensional Hamiltonians H n that 
can be factored as 



,t 



(4) 



with the periodicity condition H n+ n = H n This method 
is directly related to 2Z?-superintegrability. Indeed, if one 
of the Hamiltonians in the chain is 2D-superintegrable, it 
is possible to construct additional integrals for the whole 
family of Hamiltonians, for if [H n (x,p x ) + H / (y,p y ), Q] = 
0, then 

[H n+1 (x,p x ) + H'(y,p y ),a n+1 Qal +1 ] = 0. 

It turns out indeed that with N = 3, the first nontriv- 
ial case of dressing chains, equations |@} were solved to 
give a family of 2_D-superintegrable potential defined in 
terms of the fourth Painleve transcendant (see e.g.[l4|). 
classified as (Q.17) in 0. Since Hamiltonians satisfying 
dressing chains have interesting regularity and solvabil- 
ity properties (see 14] ), the relation between the exis- 
tence of dressing chains and integrability deserves fur- 
ther study. Here we will simply use one of the results of 
[Til to give two further examples of 2D-superintegrable 
potentials with interesting spectra. 

Let us consider the potential 



V = 



8a 4 



(x — a) 2 (x + a) 2 



with a £ M*, which is a special case of the family of 
potentials classified as (Q.17) in Q- 

This potentials behaves at infinity like the harmonic 
oscillator, but has two second-order poles at x = ±a<. 
The method of dressing chains allows us in theory to 
calculate the spectrum for this potential from that of the 
harmonic oscillator. Indeed, if we write 



= ! ' Px ~ ~2% X ~ 



b f = [p x + ^x + ihf(x) ) /V2 



(5) 



Hi 



btf 



- -l/( x - a) 


— l/(x + a) 


, we find 




■pi h 2 x 2 

fx , 


5ft 2 






2 8a 4 


a 2 






p 2 h 2 x 2 


ft 2 


ft 2 


3ft 2 


+ + 

2 8q 4 


(x — a) 2 


(x + a) 2 


4a 2 



For every eigenvector \<fr) of Hi, the function b \(j>) is ei- 
ther zero, or an eigenvector of H 2 - Conversely, if we have 
an eigenfunction \ip) of H 2 the function \ip) is zero, or 
an eigenvector of Hi. Since the spectrum and eigenvec- 
tors of Hi are well-known, we might get the impression 
that the problem for Hi is solved. This is not the case, 
since the eigenfunctions for H 2 obtained by the applica- 
tion of b on the analytic, bound states of the harmonic 
oscillator 4> n (x) have a first-order pole at x = ±a, and 
are therefore not square-integrable, unless <j) n {±a) = 0. 
Since this is the case only for <f>2, and since b<f>2 = 0, we 
can not construct a single square integrable bound state 
this way. 

Let us therefore consider a hypothetic square- 
integrable ground state ip of H 2 . We will now use an 
important result of Veselov ^4j which states that all so- 
lutions of H 2 ip = E-ip (his result applies in fact to any 
Hamiltonian of the form (Q.17)) must be meromorphic. 
We can therefore develop tp in Laurent series around 
x = a or x — —a. Since tp is square integrable and 
meromorphic, it can not diverge at any finite x, hence 
the Laurent series are in fact Taylor series. 

Let us now consider the Schrodinger equation H 2 ip = 
Eif). Since tp is nonsingular, it must be proportional to 
(x — a) 2 (x + a) 2 . Therefore b'tp, which is a eigenfunction 
of Hi, is proportional to {x — a)(x + a), and therefore 
does not diverge at x = ±a. The function wip is therefore 
an analytic eigenfunction of the harmonic oscillator. If 
b'lp — 0, we can solve directly for ip and find 



ip(x) 



exp • 



(6) 



which is not square integrable. If b^ip ^ is square inte- 
grable, it is a standard solution of the harmonic oscillator, 
and can be written as 

h(x)e~^* 2 ~ 

where h{x) is a polynomial in x. 

We have just seen that in this case b(b'ip) = Eip cannot 
be square-integrable, and therefore E = and b^tp = (p 2 , 
the second excited level of the harmonic oscillator. 

If we solve tfip — <p 2 , we find 



tp (X 



x(x 2 



a 2 )e 



-x 2 /4a 2 



2na 3 e x / 4a erf f 



This function is not square integrable because of its 
behavior at infinity and at x = ±a. 
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We have therefore proven that if ip is a square inte- 
grable eigenstate of H 2l the function b'tp cannot be a 
square integrable function. 

Since it is an analytic eigenstate of the harmonic oscil- 
lator, we know from undergraduate textbooks on quan- 
tum mechanics (i.e. |l5j ) that if b^ijj is not square inte- 
grable, it goes to infinity as a polynomial times e x / 4a . 
If that is the case, bwip = Eip cannot be square inte- 
grable at infinity unless it is zero. Since we have already 
found the two states that have E — 0, we have shown 
that the Hamiltonian H2 has no square integrable eigen- 
function. This Hamiltonian is therefore equidistant, but 
in a quite degenerate manner. Interestingly enough, if 
we consider the same Hamiltonian along the imaginary 
axis, or equivalently with ia € R, we find a real Hamilto- 
nian without singularities. We can this time find square 
integrable bound states of btf from states of the har- 
monic oscillator. We find indeed that the Hamiltonian 
bo* with a imaginary admits levels with E = 0, and with 
E = (n+3)h 2 /(2a) with n e N, but none at E = K 2 /(2a) 
or E = h 2 /a. We therefore have an example where the 
creation operator does not skip any level, but where the 
spectrum is not exactly equidistant either. 

The quantum 2D-superintegrable potentials we stud- 
ied here all have in common with the harmonic oscil- 
lator that their spectrum could be generated from a fi- 
nite number of states by the application of a creation 
operator. Their emission spectrum is highly degener- 
ate, which is exactly the kind of properties we expect 
the quantum harmonic potentials to exhibit. Moreover, 
their classical limit is the (isochronous) harmonic oscilla- 
tor. Quantum integrals yield classical integrals and it is 
therefore believed that quantum superintegrable poten- 
tials have superintegrable classical equivalents. There- 
fore quantum 2_D-superintegrable Hamiltonians are ex- 



pected to have 2D-superintegrable classical equivalents. 
This provides an additional indication that the use of 
2_D-superintegrability to generalize isochronicity or har- 
monicity is appropriate. 

All this also provides a new insight on the general 
significance of superintegrability in quantum mechanics. 
Superintegrability has already been related to supersepa- 
rability and exact solvability in quantum mechanics |16| . 
We have noticed in that quantum 2Z?-superintegrable 
potentials known today are all solutions to equations 
having the Painleve property. Moreover, the result by 
Veselov [l4|, and the relation we established between 
superintegrability and dressing chains tells us that the 
eigenstates for these Hamiltonians are often meromor- 
phic in the complex plane. Our paper shows that super- 
integrability is also related to more physical properties, 
both in classical and quantum mechanics. 

We have shown rigorously the equivalence between 
2Z)-superintegrable potentials with a minimum and 
isochronous ones. We also showed that this characteriza- 
tion of isochronous potentials could easily be transposed 
into more general contexts, such as potentials without a 
local minimum, or, more importantly, to quantum me- 
chanics. 2D-integrability can be easily generalized to 
Hamiltonians separable in different set of coordinates, or 
else in higher dimensions. Considering maximally su- 
perintegrable three-dimensional potentials of the form 
Vx(x, y) + ^2(2;), for example, allows us to deal with 2D 
isochronous potentials. The results we obtained together 
with this variety of possible generalizations demonstrates 
the fruitfulness of the approach presented here. 

The author thanks P. Wintcrnitz and C. Rousseau for 
useful discussions. This work was supported by NSERC 
of Canada. 
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